Spin Fluctuations, Interband Coupling, and Unconventional Pairing in Iron-based 

Superconductors 
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Based on an effective two-band model and using the fluctuation-exchange (FLEX) approach, we 
explore spin fluctuations and unconventional superconducting pairing in Fe-based layer superconduc- 
tors. It is elaborated that one type of interband antiferromagnetic (AF) spin fluctuation stems from 
the interband Coulomb repulsion, while the other type of intraband AF spin fluctuation originates 
from the intraband Coulomb repulsion. Due to the Fermi-surface topology, a spin-singlet extended 
s-wave superconducting state is more favorable than the nodal dxy-wave state if the interband AF 
spin fluctuation is more significant than the intraband one, otherwise vice versa. It is also revealed 
that the effective interband coupling plays an important role in the intraband pairings, which is a 
distinct feature of the present two-band system. 
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PACS numbers: 74.20.Mn, 74.20.Rp, 74.90.+n 



The recent discovery of superconductivity with higher 
transition temperatures in the family of iron-based mate- 
rials [l[ has stimulated enormous research interests both 




11 1 and theoret- 



20, 21 



experimen tally 
ically [IT ' 

In particular, the origin and nature of superconduc- 
tivity and spin density wave (SDW) ordering observed 
in these materials have been paid considerable atten- 
tion @, i, i, El El El EI EI El EI III [II . Currently 



available experimental data suggested that the supercon- 
ducting pairing state exhibit nodal behaviors [4j, while 
preliminary theoretical arguements/analyses indicated 
the pairing possibilities of an extended s-wave (either 



18], 

M 



without |14( or with nodes [l5j), a nodal d-wave IE 
a spin-triplet s-wave 17j, and a spin-triplet p-wave 
all of them are based on the scenario that spin fluctu- 
ations induce the superconductivity in this kind of sys- 
tems (with antiferromagnetic fluctuations being respon- 
sible for the former two spin-singlet pairings while ferro- 
magnetic origin for the latter two spin-triplet pairings). 
Therefore, systematic and profound theoretical investi- 
gations on spin fluctuations and their relationship with 
the superconducting pairing are significant and of current 
interest. 



This new family of superconductors has a layered 
structure, where the FeAs layer is experimentally sug- 
gested to be responsible for the superconductivit y [l| , |2|, 
IE IE H, S H • The LDA band calculations [H, EE in- 
dicate that there are five bands intersect the Fermi level 
in the folded Brillouin zone (BZ), in which four bands 
are quasi-two-dimensional. Therefore, in the representa- 
tion of the unfolded (or extended) BZ, two bands may 
be able to reproduce the main features of the four Fermi 
pockets after folding. In this paper, we employ an effec- 
tive two-band model Bamiltonian [2(| to explore the low 
energy excitation physics including spin fluctuations and 



superconducting pairing with the FLEX approach [251 ] . 
It is illustrated that one type of commensurate AF spin 
fluctuations stems from the interband Coulomb repul- 
sion associated with the nesting between the electron and 
hole Fermi pockets, while the other type of intraband AF 
spin fluctuation originates from the intraband Coulomb 
repulsion. Due to the Fermi-surface topology, a spin- 
singlet extended s-wave superconducting state is more 
favorable than the nodal dxy-wave state as the inter- 
band spin fluctuation is significant, otherwise vice versa. 
It is also elaborated that the effective interband coupling 
is enhanced by the interband AF spin fluctuation and 
plays an important role in the intraband pairings. 

We start from an effective two-band model Hamilto- 



H — Hq + Hi, 



where Hq is given by 



H 



hi a 



ei(k)cl 



klrr 



Ckla- 
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(2) 



Here Ckia denotes the band-electron annihilation operator 
with the wave vector k, spin a in the band I (I = 1, 2). In 
the present work, the two energy bands denote the hole 
band (band-1) and electron band (band-2), with their 
dispersions being approximated by £1,2 (fc) = [Cxz + £j/z T 
V^-^+^lA where £ xz (k) 
2i 2 cosfc y — 4±3 cos k x cos k y , £ yz {k) = 
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2ti cos k 



4^3 cos k x cos k y , e(k) 



— 2ti cos k x 
= —2^2 cos fc^ 
At 4 sin k x sin k 



addressed in Refs. [20j, l26|. Here, £ X z{k) and £, yz (k) may 
be understood as the iron d xz and d yz orbit-dispersions 
while e(k) as the hybridization of the two orbits. To 
produce better the topology of the Fermi surface and 
band features of the LDA calculations 13, IE 3 , we set 
t x = -1.0, t 2 = 1.5, t 3 = -1.2, t 4 = -0.95, n = 0.74, 1.7 
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FIG. 1: (Color online) (a) and (b): The Fermi surfaces of 
the two-band model in the extended Brillouin zone (1 Fe per 
cell) for /j,a — 0.74 and fib = 1.7, where the thin dashed line 
denotes the folded Brillouin zone (2 Fe per unit cell) and the 
arrow represents the nesting wave vector (see text), (c) The 
corresponding band structure (energy is in the unit of ti) in 
the folded Brillouin zone. 



(in units of |ii|), which gives rise to the electron Fermi 
pockets and the hole Fermi pockets (being respectively 
denoted by the green and red lines in Fig.l (a)) as well 
as the band structure (Fig. 1(c)). 

The interacting term Hi nt consists of the 
effective intraband 



interaction 



Hla-Hlt 



Coulomb 
cua'Cua, the effective interband 



Coulomb interaction (U'/2) J2i i^v a a > c 
the Hund's coupling 
, and the interband 
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Gil' a 

pair-hopping term 

J'Tsiw^'cliAiaiCa'a'Ca'^ where the *" site is 
defined on the reduced lattice (one Fe per cell). 

Experimental data 

@, H Hfl, ED] indicated that the 
undoped material LaOFeAs behaves like a semimetal, 
and exhibits the itinerant antifcrromagnctism. Thus it 
is reasonable to consider the Coulomb interaction to be 
intermediate in this system. In this sense, the FLEX 
approach [25| appears to be an adequate method. In 
this approach, the spin/charge fluctuations and the elec- 
tron spectra are determined self-consistently by solv- 
ing the Dyson's equation with a primary bubble- and 
ladder-type effective interaction. For the two-band sys- 
tem, the Green's function and the self-energy are ex- 
pressed as the 2x2 matrices, satisfying the Dyson equa- 
tion: G(k) = [iu n I — e(k) — S(fc)] -1 , with en = ei, 
£22 = £2 and £12 = £21 = 0. The self-energy reads 
Smn(fc) = TjY. q Y.^ v ^m, vn {q)G l _ lv {k - q), where the 
effective interaction V is a 4 x 4 matrix (with the basis 
(|11), |22), |12), |21») given by [H 



V„ 



with 



,(«) = 



[^U s x s (q)U s + lu c x c (q)U c + lu s (3) 



2 V 



+ u c )t(u s + u c )U 



X s (q) = [I - x(q)U s }~ x(q), T{q) = [I + x(q)UT l x(q) 

(4) 

as the spin and charge fluctuations. The irreducible sus- 
ceptibility is X^rn, m (q) = ~§J2k G "^i k + q)Gmn(k), 



and the interaction vertex reads 
U sl 



U s = 



u 



s2 



u c = 



u cl 
u c2 



(5) 



n and 2 J otherwise, U^ n 



where U^ n — U for m 
for m = n and 2J' otherwise 
2U' - 2J otherwise, U^ n = 
2 J' otherwise. In the above equations, k 



U£ n = U for m 
-U' + 4J for m 



U' 
n and 
n and 
(k,oj n ) and 



q — (q,iv n ) are used, with uj m the Matsubara frequency 
T the temperature, and N the lattice site number. 

The above equations form a closed set of equations 
and can be solved self-consistently to get the renormal- 
ized Green's function in the presence of the interaction 
Hi n t. After obtaining G(k), we can look into the super- 
conducting instability and the gap symmetry from the 
following Eliashberg equation [28f 



AA mn (fc) 



T 
~N 



EEE r »U^) ( 6 ) 

xG af i(k - q)Gp u {q - fc)A M ^(fc - q) 



with the pairing potential being given by T s (q) = 
%U s X s (q)U s - \U c X c {q)U c + \{U S + U c ) and f'(g) = 
-\U s X s {q)U s ~ lU c x c {q)U c + W s + U c ) for the spin- 
singlet and spin-triplet states (F'(g) is the same for the 
pairing spin projection S z — ±1,0), respectively. The 
eigenvalue A — > 1 when the supercondcuting transition 
temperature T c is reached. It is worth indicating that 
the interband Cooper pairing gap function A12 is decou- 
pled from the equation of intraband Cooper pairing gap 
functions A;; (I = 1, 2) and is vanishingly small for the 
present Fermi pockets pattern. 

Numerical calculations are carried out with 32 x 32 k- 
meshes in the extended BZ and 1024 Matsubara frequen- 
cies. The analytic continuation to the real frequency is 
carried out with the usual Pade approximant. As for the 
interaction parameters, we note that a set of parameters 
with U = 0.2 - 0.5 bandwidth and J w 0.09 bandwidth 
was used in the literature [l6j]. Here we choose 8 sets of 
representative parameters: (U, C/')=(6.5, 3.5) and (5.5, 
4.0) for J' — 1.0 &0.5 and fj, = 0.74 &1.7 with J = 1.0. 

Let us first address the static spin susceptibility 
X s {oj = 0) (x c {u = 0) < x s (to = 0), not shown 
here.). Fig. 2 presents the physical spin susceptibility 
Xp h = Y. mn X S mn,mm its intraband components x^ Xin 
and the interband one X12 m the extended BZ with U = 
6.5, U' = 3.5, J = J' = 1.0, and /x = 0.74. (For brevity 
Xmn = Xmn mn 1S used hereafter.) The physical spin sus- 
ceptibility displays two sets of dominant peaks, with one 
around (ir, 0) and its symmetric points in the extended 
BZ, and the other around (0.57T, 0.57r) and its symmetric 
points. The commensurate AF spin fluctuation around 
(tt, 0) comes from the interband Coulomb interaction as- 
sociated with the nesting between the hole and electron 
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FIG. 2: (Color online) The fc-dependence of the static spin 
susceptibility for U = 6.5, U' = 3.5, J = J' = 1, and /j, = 0.74 
at temperature T = 0.01. (a) The physical spin susceptibility 
(see text), (b)-(d) The components of the spin susceptibility 
X22, Xii, an d Xl2, respectively. 



pockets [18| (Fig.l (a)), which is clearly seen from the in- 
terband component \i2 shown in Fig. 2(d). The appear- 
ance of this type of AF spin fluctuation around (±7r, 0) 
and (0, ±7r) is in goo d ag reement with the neutron scat- 
tering measurements [H 



and here we refer to it as 
the interband AF spin fluctuation. More intriguingly, 
the other type of intraband spin fluctuation is seen to 
peak around (0.57T, 0.57r) in the components xfi and X22J 
which is mainly induced by the intraband Coulomb inter- 
action U associated with the approximate nesting prop- 
erty within the renormalized Fermi pocket around the 
(ir, 7r) point (not shown here). Notably, the peak posi- 
tion corresponds to [0,7r] (and [ir,0]) in the folded BZ 
and thus implies the emergence of a new component of 
"stripe" -type AF spin fluctuation in the primary lattice 
(with 2 Fe ions per unit cell) , which could be referred to 
as the intraband spin fluctuation and may be detected di- 
rectly by future neutron scattering experiments on single 
crystal samples. 

The most favorable superconducting pairing symme- 
try at a fixed temperature is determined by solving the 
Eliashberg equation with the maximum eigenvalue. The 
calculated maximum eigenvalues (for various possible 
pairing symmetries) versus temperature are plotted in 
Fig. 3 (a). Firstly, one can see that the eigenvalue for the 
spin-triplet p-wave state is much smaller than those of 
the spin-singlet state, and at the mean time exhibits a 
flat temperature dependence. Therefore, we can safely 
rule out the possibility of the spin-triplet state in the 
present model calculation. In the spin-singlet channel, 
the eigenvalue of dxF-wave state is larger than that of 
the s-wave state, and in particular the former increases 
rather rapidly with decreasing temperature. In view of 
this tendency, although the maximum eigenvalue A = 1 
has not been reached yet, it is reasonable to consider the 
spin-singlet d^y-wave to be the most favorable state in 
this set of parameters. The calculated fc-space structure 



FIG. 3: (Color online) (a): Temperature dependence of the 
maximum eigenvalues, (b) and (c): fc-dependence of the gap 
functions A2,i(fc) corresponding to the largest eigenvalue, for 
the same set of parameters as those in Fig. 2 at temperature 
T = 0.01. The solid diamond is the folded BZ and the dashed 
circle denotes the Fermi pocket schematically. 



of the gap functions for both the electron and hole bands 
in the extended BZ are depicted in Figs. 3(b) and (c), 
respectively. It is seen that the pairing symmetries in 
both bands are of dxy-wave, namely, the gap function 
A«(k) « A^7; k with j tk 2 sin k x sin ky , where (k x , 
ky) is the wave vector denoted in the folded BZ. Inter- 
estingly, we note that the gap magnitude in the electron- 
band is significantly larger than that in the hole-band. 
This feature could also be understood as follows. Eq.(6) 
may approximately be rewritten as 



(7) 



for u) = 0, where Kw = ^~Jk k' ~ k') with the ef- 

fective intraband pairing potential and interband cou- 
plingas IV = ~\G vv \ 2 [(U 2 + AJ 2 )x s w 5w + 4tf'J'xJ, (1 - 
*H')]7/k7/'k'/(-^Z)k7fk)- Since ^22 > Kn (mainly due 
to the result X22 ^ Xii); the gap amplitude of electron 
band (band 2) is mainly determined from the intraband 
pairing of itself, while that of the hole-band depends 
mainly on the interband coupling from the electron band. 
From the expression of Vu> , it is elucidated that the intra- 
band spin AF fluctuation leads to the intraband pairings, 
while the interband AF spin fluctuation enhances the ef- 
fective interband coupling (associated with the interband 
coupling factor U'J') and thus the intraband pairings of 
both bands [29(. The above conclusions are unchanged 
for /i = 1.7 or J' — 0.5, where the intraband spin fluctu- 
ation has an even more impact on the pairing. However, 
if the interband AF spin fluctuation is dominant over the 
intraband one, as seen in Fig. 4(a) for U = 5.5, U' = 4.0, 
,/ = J' = 1, and /x = 0.74, an extended s-wave state 
(Figs. 4 (b) and (c)) would be more favorable than the 
d-wave one. A similar conclusion has also be obtained in 
a recent renormaliztion group study [3(| • 

Because the spin fluctuation x' s is stronger than the 
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charge fluctuation x° ( n ot shown here), the pairing inter- 
action in the spin-singlet channel is positive {T s (q) > 0). 
Consequently, once the gap function satisfies the condi- 
tion A ; (k)A;/(k + Q) < 0, where Q is the wave-vector 
around which the pairing interaction peaks, we are able 
to obtain the largest eigenvalue solution of the Eliash- 
berg equation. Focusing merely on the interband AF 
spin fluctuation, the above condition leads to two candi- 
dates of pairing symmetry on the singlet channel. One 
is the extended s-wave, with the gap function of each 
Fermi pocket having the same sign, while changing the 
sign between the electron and hole pockets, as shown in 
Fig. 4(b) and (c). The other is the <i-wave as shown in 
Fig. 3(b) and (c). Actually, which one is more favored 
depends mainly on the existence of the intraband spin 
fluctuation caused by the Coulomb interaction U. For a 
larger U and the approximate nesting within the renor- 
malized Fermi pocket, this spin fluctuation that peaks at 
(0.57T, 0.57r) emerges (see Fig. 2(b) and (c)), which leads 
to the gap function to change its sign within each Fermi 
pocket and thus induces a nodal dxy-wave pairing. In 
contrast, if the intraband Coulomb interaction U is rel- 
atively weak, such that the intraband spin fluctuation is 
not significant, the extended s-wave state would energet- 
ically be favored as it opens a full gap around the Fermi 
pockets [3l| . We think that this kind of connection be- 
tween the peak structure of the spin response and the 
pairing symmetry established here is useful for probing 
the superconducting pairing symmetry by measuring the 
fc-dependence of spin fluctuations in neutron scattering. 

Finally, we note that the LDA energy bands in this 
system exhibits a more complex structure 12, 14, [l5j |. 
In the present two-band model, the d xz and d yz orbitals 
are chosen as these two orbitals have the lar gest weights 
to the energy bands crossing the Fermi level [32|. While, 
the d xy orbital contributes also a weight to the energy 
band along the T — M direction [32| , and particularly in 
the unfolded Brillouin zone one of the hole Fermi pockets 
is displaced from the (0, 0) point to the (tt, tt) point as 
depicted in Figs. 1(a) and (b), which is actually different 
from that in a more realistic energy band structure [HI]. 
However, we wish to indicate that these two deficiencies 
do not affect meaningfully the nesting properties of the 
Fermi pockets, and thus the results/conclusions obtained 
above are still unchanged, at least qualitatively. 

In summary, based on an effective two-band model and 
using the fluctuation-exchange approach, we have inves- 
tigated spin fluctuations and superconductivity as well 
as the interband coupling in iron-based layered super- 
conductors. We have elaborated that one type of com- 
mensurate AF spin fluctuation comes from the interband 
Coulomb interaction associated with the nesting between 
the hole and electron Fermi pockets, while the other type 
of intraband AF spin fluctuation originates from the in- 
traband Coulomb repulsion. We have elucidated that, if 
the interband AF spin fluctuation plays a significant role, 




FIG. 4: (Color online) (a): The physical spin susceptibility, 
(b) and (c): fc-dependence of the gap functions A2,i(fe) corre- 
sponding to the largest eigenvalue, for U = 5.5, U' = 4.0, J = 
J' = 1.0, and fi = 0.74 at temperature T = 0.01. 



this fluctuation leads to the pairing with the spin-singlet 
extended s-wave being the most favorable state. Other- 
wise, the pairing is mainly determined by the intraband 
AF spin fluctuation, with dxv-wave symmetry. 

Note added-(\) Recently, we note that the angle- 
resolved photoemission spectroscopy measurement on the 
(Ba,K)Fe 2 As 2 superconductor [33J] indicated the pres- 
ence of nodeless gaps, which might imply the interband 
spin fluctuation plays likely a dominant role at least in 
these iron-based superconductors, (ii) After the work 
was posted on the e-Print archive, we note a similar work 
done by X.-L. Qi et al in [34j |. 
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